Quantum Walk of Two Interacting Bosons 
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We study the effect of interactions on the bosonic two-particle quantum walk and its corresponding 
spatial correlations. The combined effect of interactions and Hanbury-Brown Twiss interference 
results in unique spatial correlations which depend on the strength of the interaction, but not on its 
sign. The results are explained in light of the two-particle spectrum and the physics of attractively 
and repulsively bound pairs. We experimentally measure the weak interaction limit of these effects 
in nonlinear photonic lattices. Finally, we discuss an experimental approach to observe the strong 
interaction limit using single atoms in optical lattices. 
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Introduction.- Understanding highly correlated many 
body systems remains both an experimental and theoreti- 
cal challenge. While there is a rather good understanding 
of weakly interacting systems, problems involving strong 
interactions are in general harder to address. 

Recently a new approach to the study of quantum 
dynamics became experimentally accessible through the 
study of Quantum Walks (QWs) in lattice potentials [lj]. 
Quantum walks are the quantum counterparts of classi- 
cal random walks on discrete lattice: A quantum particle 
is initially placed at a particular site of a lattice and then 
tunnels to neighboring sites with equal probability am- 
plitude. This basic "step" is repeated, but in contrast to 
the classical case quantum mechanical interference leads 
to distinctively different dynamics. For example, in pe- 
riodic lattices the wavefunction width grows ballistically, 
while in the classical case the expansion is diffusive. 

QWs receive increasing attention due to their relation 
to various physical and bio- physical processes [H-Qj and 
their possible use as a primitive for quantum computation 
algorithms Theoretically, QWs were studied for the 
single particle case [lj]. Initial experiments studied the 
physics of single particles by using either classical waves 
[H, single photons 0, 0], or single atoms [H, Mov- 
ing from one to two non-interacting particles it has re- 
cently been shown that indistinguishable quantum walk- 
ers can develop non-trivial correlations due to Hanbury 
Brown- Twiss (HBT) interferences 

[IMS- 

Yet, very lit- 
tle is known on the effect of interactions on the dynamics 
of the few-body QW [13| . As new systems emerge that 
can accommodate such experiments [1, 0], a systematic 
study of this problem starting at small particle numbers 
may offer a "bottom up" approach in the general thrive 
to understand dynamical quantum many-body systems. 

In this letter we study the effect of inter-particle inter- 
actions on the two-particle quantum walk and the result- 
ing spatial correlation. We consider two bosons, each ini- 
tially localized on a single lattice site, undergoing a QW 
simultaneously. We find that the interplay between in- 
teractions and quantum two-particle (HBT) interference 
gives rise to fermion-like spatial correlations between the 
particles. Interestingly, the correlations depend on the 
strength of the interaction but not on whether it is at- 




Figure 1: (color online). Two-particle Quantum Walk, (a) 
An illustration of two identical atoms initially placed on two 
sites of an optical lattice, and allowed to tunnel to neighboring 
sites and interfere, (b) The evolution of the particle-density 
for the case of a single atom quantum walk, (c) The evolution 
of the density for a two-atom quantum walk. This density is 
only weakly affected by interactions, while the two particle 
correlations are strongly modified. 



tractive or repulsive. We explain the observed correla- 
tions by calculating the two-particle spectrum, and in- 
terpret our results in light of the physics of attractively 
and repulsively bound pairs (13 . fla ] . We then present an 
experimental observation of the weak interaction limit of 
these effects in nonlinear photonic lattices, and outline an 
experimental approach to observe the strong interaction 
limit using single atoms in optical lattices. 

QW of two interacting particles.- We start by calculat- 
ing the QW of two interacting particles. We consider the 
one dimensional Bose-Hubbard model: 
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where aj^ (a rn ) is the creation (annihilation) operator for 



a particle at site m, 



is the corresponding 



number operator, J is the tunneling amplitude between 
nearest neighbors and U is the on-site interaction energy 
which can be attractive (negative) or repulsive (positive) . 

We study the QW of two indistinguishable particles, 
each initially localized on a single site in a periodic lattice. 
We consider two different initial conditions: One in which 
the two particles are localized at adjacent lattice sites 
\ipinitiai) = o|ag|0), and a second in which the particles 



are initially placed at the same site, \rp. 
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Our focus lies on the particle-density n r (t)=(a\.a r \ and 
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Figure 2: (color online). Two-particle correlations for in- 
teracting quantum walkers. Left column: correlations after 
propagation time T = 4, where initially the two particles are 
placed at the same site, \ipimUal) = ( a o) 2 |0). (a) For zero in- 
teractions, the two-particle correlation shows no interference 
[ic| |. (b)-(f) As the interaction is increased the correlations 
shows the formation of bound pairs, while the density distri- 
bution (shown on the right of each plot) becomes increasingly 
localized. Right column: similar results for an initial condi- 
tion in which the two particles are placed at adjacent sites 
\ipinitial) = a ! a o|0). (g) At U — the correlations shows spa- 
tial bunching, (h)-(j) The correlations change as the interac- 
tion \U\ is increased, while the density is only weakly affected, 
(k) At strong interactions the correlation is transformed to 
spatial anti-bunching, similar to the correlation that would 
be exhibited by two non-interacting fermions initially placed 
in the same configuration (1). All results are identical for both 
attractive and repulsive interactions. 



on the two-particle correlation T q , r (t) = ^aJa|a r o g ^ 
which are calculated after an evolution time T for differ- 
ent values of the interaction U. At all stages the particles 
are far from the lattice boundaries. 

The results for U — correspond to the results re- 
ported in (ToL 11 1 : when the two bosons start the QW 



at the same site, after propagation each particle can be 
found on either side of the site of origin, reflected in the 
four symmetric peaks in the correlation matrix inset for 
Fig. 2a. When the particles are initially placed at adja- 
cent sites, HBT interference results in spatial bunching, 
and the two particles propagate together either to the left 
side of the distribution (peak at the bottom left corner 
of the correlation matrix in [5^) or to the right (top right 
corner). Other initial conditions, in which the particles 
are further separated in space result in more complicated 
correlation patterns [Icj . Such initial state and the re- 
sults of interactions will be discussed in (lfj| . 

Let us now turn to the discussion of interaction effects. 
Fig. [5] b-e show the results for increasing repulsive inter- 
action U for in the case of two bosons initially localized at 
the same site. The spatial correlations show that as \U\ 
increases the two particles tend to propagate as a pair, 
while the density distribution becomes localized. Fig. [5] 
g-1 shows the results for the case in which the particles 
are initially places at different sites ipiniUai) = ai&o|0). 
Here, the particle density depends only weakly on U, 
but the two-particle correlation undergoes a fundamen- 
tal change: the spatial bunching effects which occur in 
the non-interacting case gradually transforms to spatial 
anti-bunching (Fig. [5] k). For large values of the in- 
teraction strength \U\, the correlation between the two 
bosons becomes very similar to the correlation exhibited 
by two non-interacting fermions, prepared in the same 
initial configuration (compare Fig. [2] k and 1). The non- 
interacting fermionic and the interacting bosonic matri- 
ces become identical at the limit of \U\ — > oo, while the 
density becomes identical to the one for U = 0. An inter- 
esting result is that in both cases the effect of interactions 
does not depend on the sign of U; it is identical for both 
attractive and repulsive interactions. We note that for 
initial conditions in which the two particles are further 
separated in space interactions also drive the system to- 
wards fermion-like correlations, only that now they have 
a more complicated spatial structure- see [l6j for addi- 
tional experimental and theoretical results. 

The two particle spectrum.- To understand these re- 
sults we consider the two-particle spectrum of the sys- 
tem, as shown in Fig. [3] for two particles on a lattice with 
AI = 29 sites. Each of the two-particle eigenfunctions can 
be written as ^(ri, r-z) where n and r-i are the positions 
of the two particles on the lattice. Introducing the center 
of mass coordinate, R = (ri+r2)/2 and the relative coor- 
dinate, r = n — ra, we can solve the Schrodinger equation 
with the ansatz ^(ri,r2) = exp(iKR)tpK{'r), where K is 
the quasi-momentum of the center of mass motion and 
iPk(t) is the pair wavefunction (l5| . 

For finite interaction strength U the spectrum sepa- 
rates into two bands. The main part of the spectrum, 
containing [M(M— l)]/2 eigenstates, consists of scatter- 
ing states having low probability at r = 0, whose en- 
ergy is given by the non-interacting part of the Hamil- 
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tonian. The smaller part of the spectrum (M eigen- 
states) consists of states ip^(r) which have a large prob- 
ability for two particles to occupy the same site, i.e., 

\^k(°)\ 2 1 ( u °°) ( see insets in Fig. EH). This 
mini-band has higher or lower energies than the main 
part of the spectrum, depending on the sign of the inter- 
action; nevertheless, the spatial probability distribution 
of the two-particle eigenstates is identical. 

Using this picture, it is possible to explain the results 
in Fig. [21 An initial state in which the two particles 
occupy the same site with strong attractive or repulsive 
interaction, will mostly contain two-particle states from 
the smaller mini- band. As a result, the two particles 
will remain bound as described by Winkler et al. in [14 1 
(see Fig. [2] g-1). A complementary process happens if 
the particles initially occupy different sites. This initial 
condition excites mainly scattering states from the main 
part of the spectrum. As a result, the particles have low 
probability to be found at the same site throughout the 
evolution, and will not show bunching. 

Let us now turn to the case of strong interactions 
\U\ 3> J. Our goal is to understand the "fermioniza- 
tion" as observed in the correlator T g:r (t) for an initial 
state in which the particles are found at different sites. 
We start by noting that by focusing on the scattering 
states we can describe the Hamiltonian (J2| using hard- 
core bosons, where doubly occupied sites are eliminated 
from the Hilbert space. Formally we replace the bosonic 
with spin-1/2 operators: a] n — > S+, a m — > S 1 " . 

Next, we use a standard mapping from spin-1/2 to 
fermionic operators f m , Let us review the essen- 

tial steps of this mapping to understand the "fermionic" 
behavior of T qtT .(t). Spin-1/2 and fermionic operators 
share the local property (/tj 2 = = (5+) 2 = {S m ) 2 = 
0. However, spins on different sites commute, whereas 
fermions pick up a minus sign. In the sought mapping one 
corrects for this via the Jordan- Wigner string exp(i</> m ): 
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It is now straight-forward to check that for T qtr (t) the 
Jordan- Wigner string drops out. Hence, the correlation 
for hard-core bosons are identical to the ones obtained 
for non-interacting fermions in accordance with our ob- 
servation in Fig. [2] 

Experimental results.- The case of the two-photon 
quantum walk and the resulting HBT correlations were 
considered in [lj| and observed in [ll[, in a system of 
waveguide lattices. This system is described by an equa- 
tion identical to Eq. 1, only that for single photons in- 
teractions are negligible, i.e. U = 0. Ref. [10( presented 
also a measurement of correlation for classical, thermal 
light waves, analogous to the intensity correlations pre- 
dicted by the original HBT work 12j. That experiment 
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Figure 3: (color online). Spectrum and two-particle eigen- 
modes of Eq. [2]for U = 8. The spectrum is separated into two 
mini-bands. The higher band consists of bound pair states in 
which the two particles only occupy the same sites (non-zero 
values mainly on the diagonal r = q, top inset), while the 
lower band consists of states in which the two particles have 
small probability to occupy the same site (lower inset). The 
gap is proportional to the interaction strength U. Attractive 
interactions U < 0, yield an identical yet inverted spectrum. 



showed that the results for thermal inputs captures some 
(but not all) aspects of the correlations predicted for the 
quantum, two-particle case. 

The waveguide lattice used in that experiment, de- 
scribed in detail in 0, [To| , has been shown to supp ort 
nonlinear effects for high intensity classical light ll8|. In 
this case the system is described by the Hamiltonian: 



(l.m) m 



(2) 



which is the classical or mean-field limit of Eq. 1. It is 
therefore reasonable to presume that the classical corre- 
lations for nonlinear thermal waves in this system will 
correspond to the results for two interacting quantum 
particles. Indeed, In the experiments described below we 
find that in the limit of weak interactions the measured 
classical correlation for nonlinear thermal waves are sim- 
ilar to the predicted quantum correlations. 

In Fig. [4] we present experimental measurements for 
intensity correlations obtained using nonlinear thermal 
waves in \ipinitial) = a I a o|0)- Detailed numerical results 
are presented in Fig [5] In all figures we compare the cor- 
relation fluctuations: F qr (t) = (aJaJa r o g ) — \ (o-\o-q) • 



(ala r ) 



\n q ■ n r which are a better basis for 



comparison between the quantum and classical (thermal) 
case. As the results show, for weak interactions the classi- 
cal HBT correlations follow the quantum predictions - see 
Fig. 0] and Fig. [5] a-f. However, as interactions become 
stronger the two systems diverge - while the quantum 
system exhibits a switch to fermion-like correlations, the 
classical system cannot follow, and remains with modi- 
fied, localized correlations - Fig. [5]i-l. 

Proposed cold atoms experiment.- Experimentally, the 
strong interaction limit of effects discussed above can be 
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Figure 4: (color online). Experimental measurements of 
the fluctuation in intensity correlations for nonlinear thermal 
light, (a) The fluctuations in the intensity correlations in the 
linear case, when two thermal beams are injected into two ad- 
jacent sites, corresponding to \ipi n iuai) — a i a o|0). The results 
show spatial bunching [10] ]. (b) the predictions of the quan- 
tum theory for two, non-interacting particles initially placed 
at the same locations, (c) Experimental results for nonlin- 
ear thermal waves (d) the predictions of the quantum theory 
for two interacting particles. Note the similarity between the 
classical results and the quantum prediction in both cases. 
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control interactions, for example via tuning a Feshbach 
resonance, or controlling the ratio U /J. For the two par- 
ticle case and at zero interactions, the results should cor- 
respond to those presented in [nj and observed in fill ] 
using photon pairs. However, when interactions will be 
introduced we predict the results presented above: if the 
particles are places one on top of each other they will 
tunnel as a pair and the density will become local- 
ized. If they are placed at different, not too distant sites, 
the density will show only minor changes as a function 
of U, but the two-particle correlation will change signif- 
icantly, reaching a fermion-like correlations at the limit 
of strong interactions. In the same spirit, this system 
can be used to directly measure the dynamic properties 
and correlations for large number of particles, a problem 
which quickly becomes impossible to compute. 

Conclusions.- In this letter we have considered the 
quantum dynamics of two bosons on a lattice, each ini- 
tially confined to a single site. Such dynamics with two 
or more particles can be experimentally explored in sys- 
tems such as described in ref. As the number of par- 
ticles increases the problem will become uncomputable, 
but may remain experimentally accessible. 

Note.- During the final completion of this manuscript 
we became aware of related work [20fl. 



Figure 5: (color online). Simulations comparing the fluctua- 
tions for the quantum (top panel) and the classical (bottom) 
cases, with increasing \U\ or nonlinear coefficient |7|, corre- 
spondingly, for \ipi„itiai) = a i a o|0). The two sets of results 
look similar for weak interactions (a-c). However, beyond 
\U\ — 1.5 (d-g) the two results diverge - the quantum cor- 
relations transform to fermionic-like correlations, while the 
classical correlations become increasingly localized. 

observed using techniques that recently became available. 
The experimental requirements include the ability to: (i) 
initially localize exactly two quantum particles at two 
pre-determined lattice sites, (ii) to allow these particles 
to freely tunnel and exhibit a QW, (iii) to control the in- 
teraction strength (or the interaction-to-tunneling ratio), 
and (iv) to image single particles in the lattice sites after 
some evolution time. An example of a system that can 
accommodate such experiments was recently presented 
in In this system, the authors placed single atoms at 
selected sites and allowed them to tunnel, exhibiting, on 
the ensemble average, the dynamics of continuous time 
QWs [H,i|. Using a similar approach, we propose start- 
ing with an ensemble of two (and in principle N) atoms 
separated by several sites. The density after time T can 
be measured in the same manner as in and the two 
(or N) particle correlation can be directly assessed from 
the raw data. An important aspect will be the ability to 
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